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Abstract 

Background: Mechanical system with imperfect or misalignment of X-ray phase 
contrast imaging (XPCI) components causes projection data misplaced, and thus 
result in the reconstructed slice images of computed tomography (CT) blurred or 
with edge artifacts. So the features of biological microstructures to be investigated 
are destroyed unexpectedly, and the spatial resolution of XPCI image is decreased. 
It makes data correction an essential pre-processing step for CT reconstruction of XPCI. 

Methods: To remove unexpected blurs and edge artifacts, a mathematics model for 
in-line XPCI is built by considering primary geometric parameters which include a 
rotation angle and a shift variant in this paper. Optimal geometric parameters are 
achieved by finding the solution of a maximization problem. And an iterative approach 
is employed to solve the maximization problem by using a two-step scheme which 
includes performing a composite geometric transformation and then following a linear 
regression process. After applying the geometric transformation with optimal 
parameters to projection data, standard filtered back-projection algorithm is used to 
reconstruct CT slice images. 

Results: Numerical experiments were carried out on both synthetic and real in-line 
XPCI datasets. Experimental results demonstrate that the proposed method improves 
CT image quality by removing both blurring and edge artifacts at the same time 
compared to existing correction methods. 

Conclusions: The method proposed in this paper provides an effective projection data 
correction scheme and significantly improves the image quality by removing both 
blurring and edge artifacts at the same time for in-line XPCI. It is easy to implement and 
can also be extended to other XPCI techniques. 

Keywords: Projection data correction, Tomography reconstruction, In-line X-ray phase 
contrast imaging 



Background 

In recent years, X-ray phase contrast imaging (XPCI) has attracted much attention in 
many areas, such as bio-medical imaging [1], material science [2] and paleontology [3]. 
Different from conventional absorption-based X-ray imaging, XPCI uses the informa- 
tion of wave front changes while X-ray radiation passes through a sample. It produces 
refraction-sensitive images and reveals the internal anatomy of soft tissues with diffe- 
rences in the refractive indices, which can be described in the complex form as follows 
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n = l- S + iji 



(1) 



where 6 is the decrement of the real part of the refractive index, and [3 presents the ab- 
sorption index of local attenuation. It is proved that 5 is three orders of magnitude lar- 
ger than p within the diagnostic X-ray range for human tissue [4], Thus XPCI can 
greatly enhance the contrast for soft tissues in biological sample with high spatial and 
temporal resolutions than conventional absorption-based X-ray imaging technique. 

To reveal the morphology of a thick biological sample by using XPCI, 2D projection 
images are collected by titling the sample around a fixed axis, which is perpendicular 
to the beam direction, and computed tomography (CT) technique is used to obtain 3D 
visualization of the internal structure of the sample. Usually standard filtered back- 
projection algorithm (FBP) or algebraic reconstruction technique (ART) [5] are 
employed to build CT slice images from projections. However, collected projections 
cannot be used for CT reconstruction directly without pre-processing in practice. Due 
to mechanical system with imperfect or misalignment of X-ray source, rotary stagy 
and CCD detector, an unavoidable problem encountered before CT reconstruction is 
how to align 2D projections and make the rotation axis coincident with the center line 
of each projection. Similar problem has been investigated in electron tomography 
(ET). A conventional scheme for this issue in ET is to employ gold particles or image 
derived markers to correctly align projections by computing geometric parameters 
from tracking and measuring projected positions of them [6,7]. However, this method 
needs additional components to generate accurate and measurable markers in projec- 
tion images and thus makes the imaging procedure and post processing complex. 
Some marker-free methods are also proposed to handle this issue in ET by using the 
cross-correlation (CC) of successive projection images [8] or taking specimen features 
as built-in markers for deriving geometric parameters [9,10]. These methods are based 
on pattern matching and have the drawback of accumulating alignment errors by com- 
paring successive projection images. Recently, a mass center (MC) based method 
[11,12] is proposed to align 2D projection images before reconstructing 3D tomog- 
raphy image in coherent diffraction imaging (CDI) [13,14]. This method does not use 
local features but mass center as the invisible marker. Due to only one marker can be 
used, MC-based method works well for translational alignment, but poorly for rota- 
tional alignment. And it is sensitive to background noise and thus needs a well- 
designed background subtraction scheme. Some existing software suites developed for 
XPCI reconstruction usually provide a manual and automatic integrated approach to 
handle the data correction issue before CT reconstruction. PITRE [15] employs a MC- 
based sinogram correction scheme to calibrate the location of rotation axis, and the 
user can refine the result manually. DEIRecontructor [16] calculates geometric param- 
eters by using user-selected markers. In this paper, a new automatic method is pre- 
sented to solve the data correction problem in in-line XPCI, which is a propagation- 
based phase contrast imaging technique. This method employs an iterative approach 
to properly determine optimal geometric parameters. And a two-step method, which 
includes performing a composite geometric transformation and then following a linear 
regression process, is used to find out the solution at each iteration. Numerical experi- 
ments on synthetic and real image datasets demonstrate that the proposed method 
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provides a fast and reliable image pre-processing scheme for CT reconstruction in in- 
line XPCI. 

The rest parts of this article are organized as follows: Methods section describes 
the projection data misplacement problem for CT reconstruction in in-line XPCI 
and its mathematical model with primary geometric parameters at first, and then in- 
troduces the proposed method in details. Results section shows numerical experi- 
mental results from synthetic and real datasets which are pre-processed by the 
proposed method, and followed by the study conclusion in Consclusions section. 

Methods 

Projection data correction in in-line XPCI 

XPCI uses a spatially coherent beam and an X-ray sensitive detector to acquire pro- 
jection images from a sample. To obtain 3D structural information of the sample, a 
rotary stage with high precision is used to put the sample on and a series of projec- 
tions collected while titling the sample around a rotation axis with respect to the 
base of stage. In-line XPCI [17] is a propagation-based phase contrast imaging tech- 
nique and the most simple and straight-forward method to achieve X-ray phase con- 
trast images. A synchrotron X-ray source, the sample and an X-ray sensitive detector 
are arranged in line as in conventional radiography but the detector is placed behind 
the sample in a certain distance. Figure 1 shows the schematic set-up of in-line 
XPCI. 

While reconstructing slice images from acquired projections in in-line XPCI, the 
center line of any projection is considered to be identical to the rotation axis by de- 
fault. Following this assumption, reconstruction algorithms, such as FBP, can be 
employed to build slice images correctly. However, it is usually very hard and time- 
consuming to align the components of imaging system with micrometer precision. If 
no data pre-processing step is adapted, the final reconstructed slice images will be 
blurred and with edge artifacts. 

Suppose that Y is along the center line of CCD projection and Z is along the X-ray 
beam direction. Considering CCD detector and X-ray beam source are fixed during 
the whole imaging process, we use a model with four geometric parameters to de- 
scribe the space relationship between the center line of projection and the rotation 
axis at tilt angle a as illustrated in Figure 2. Here 6 is the in-plane flip angle which 
specifies the slope of the rotation axis related to Y, and (S x , S y , 8 Z ) are offsets along X, 
Y and Z directions respectively. S y and S z are caused by the movement biases of ro- 
tary stage, and 9 and S x are mainly due to the misalignment of CCD detector and ro- 
tary stage. The movement biases of rotary stage are controlled by the precision of 
mechanical design, and thus S y and S z are usually small and can be ignored. So 6 and 
5 X are two main factors to be considered and the model is simplified as illustrated in 
Figure 3. There S x is denoted as S. The existence of two misalignment parameters re- 
sults in geometric transforms of the projection images. The parameter 6 is related to 
the rotation transform of projection images while S to translational transform along 
X. So the misalignment problem can be solved by applying inverse geometric trans- 
form to the projection images if we know the values of these geometric parameters. 

Figure 4 compares central slices reconstructed from projections of a modified 3D 
Shepp-Logan phantom [18] without and with geometric transforms by using FBP. 
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Figure 1 The schematic set-up of in-line XPCI. 



Translational transform with S = 2 pixels and rotation transform with 6=5° are 
applied to projections respectively. Compared to the result without geometric trans- 
form, translational transform makes the slice image blurred while rotational trans- 
form introduces obviously edge artifacts. 



Two-step iterative correction method 

In-line XPCI has no limitation on the range of tilt angle. Usually two projections 
with special tilt angle of 0 and 180 degree are included while collecting projections 
in in-line XPCI. Denote P 0 and P„ the two projections with special tilt angles re- 
spectively. According to the physical principle of in-line XPCI, P„ is the reflection of 
P 0 by taking the rotation axis as the mirror line. So P n should be identical to the 
image which is achieved by applying the following composite geometric transfor- 
mation to Po in homogeneous coordinates 




Wu ef al. BioMedical Engineering OnLine 2014, 13:105 
http://www.biomedical-engineering-online.eom/content/13/1/105 



Page 5 of 1 1 



Y ' 


1 i 




jL / 

i / 

/ \ 1 ■ 




" ) 
/I / 1 

/ 1 / 

1 1 / 1 

i — 1 


i , 
i 

i 

i 

i 

i 


1 
1 

1 
1 
1 



Figure 3 Geometric parameters in simplified model. 6 denotes the in-plane flip angle and 5 is the offset 
from rotation axis to the center of projection image along X. Dashed box presents the projection window. 
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where (x,y) is a point in P 0 and (x',y') is its corresponding position in P n . If the rota- 
tion angle 8 is small enough, then (x' ,y') is approximated by 



/ «y + 2 (tf-*)0 



Thus we obtain 

Po(x,y) = Pn{2{6 - By) - x,y + 2{S - x)d) 



(2) 



(3) 



The above equation shows that the optimal geometric parameters can be achieved by 
solving the following minimization problem when P 0 and P„ are available 



© © (D 

(b) (c) (d) 



Figure 4 Central slice reconstructed from projections of modified 3D Shepp-Logan phantom. From 
left to right: (a) the original, (b) reconstructed without geometric transform, (c) reconstructed with 6 = 2 
pixels and (d) reconstructed with 0 = -5°. 
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(0*,<T) = arg min jj{P 0 {x,y) - P n {2(8 - 8y) - x,y + 2{S - x)8)) 2 dxdy (4) 

8,8 

Rewrite Eq. (4) by discarding the constant terms as follows 
(6>*, 5*) = arg max fl{P 0 {x,y)P n (2(S - 8y)-x,y+2(S - x)8))dxdy (5) 

To find the solution of the above optimization problem, we employ an iterative 
scheme. At the Ath iteration, the basic step is given by the following iteration 

(d k ,8 k ^ = arg max j (p 0 (x,y)P n (2(8 - By) - x,y + 2^8 k .i - x^dk-ij^dxdy (6) 

The optimizer of Eq. (6) is achieved by using a two-step method which includes per- 
forming a composite geometric transformation 1 and then following a linear regres- 
sion process. 

In the Ath step, 7* ~ 1 is applied to P„ to rotate it with the angle - 28 k _ i and then re- 
flect it. Here T* ~ 1 is defined as follows 
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P k n is obtained by applying T* 1 to P n . It can be proved that we have 

Pl(x,y) = P„(2(8 -6y)-x,y+ 2{8-x)6) (7) 
Then we seek to find the solution of the following optimization problem 

(e k ~8 k ) = arg maxjr(p 0 (*,y)l£(2((* - 5 k -i) - (o - Ok-i)y) +x iy ))dxdy (8) 

Considering the cross-correlation of P 0 and P k n along X 

(P 0 * (t) = \Po(x,y)P k n (t + x,y)dx, (9) 
Eq. (8) can be rewritten as follows 

(e k ,Sk) = arg max/(P 0 * i*) (2^8 - S k -i) - (b - d k -i)y))dy (10) 

The above equation shows that the optimizer of Eq. (8) can be achieved by comput- 
ing the cross-correlations of P 0 and P* row by row, and fitting the positions t m (y), 
which are corresponding to the maximum of correlation coefficients, with the following 
linear model 



t m {y) = 2(8 - 8 k -i) - 2(e - 8 k .^y 



(11) 



The stop criterion of the iterative method is set as both 8 and r5 are less than speci- 
fied tolerance values respectively. T° is set to the identity matrix. 
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Results 

To evaluate the proposed method, we carried out experiments on both synthetic 
and real in-line XPCI data. MC-based sinogram correction scheme and the proposed 
method have been implemented using MATLAB (The MathWorks, Inc., Natick, 
MA, USA). Inverse Radon transform with Ram-Lak filter was taken as CT recons- 
truction algorithm applied to corrected data. Data processing for this paper was car- 
ried out on a Dell workstation system with a 2.4GHz Intel Core i5 processor and 
8GB memory. 

Simulation on synthetic data 

For the synthetic case, we used projections of modified 3D Shepp-Long (S-L) phan- 
tom in Figure 4 of size 256x256x256 as a test dataset. These projections were ob- 
tained by using Radon transform after applying geometric transformation with 
specified parameters to the 3D S-L phantom. 2 pixels shift along axis X was applied 
to the rotation axis while flip angle was set to -5 degree. Tilt angle ranged from 0° to 
180° and 181 projections were built in the simulation with one projection per degree. 
Figure 5 shows the results of 4 middle CT slices reconstructed without correction 
(NC), with MC-based sinogram correction (MC), and with proposed geometric cor- 
rection method(GC) respectively. By comparing these slice images, we confirm that 
MC-based sinogram correction scheme reduce the blur effect but cannot remove 

r \ 



©O©0 
©©©© 
0©©© 
0©©© 



Figure 5 Four CT slice images of modified 3D Shepp-Logan phantom reconstructed by FBP. From 
top to bottom are the results of references, without correction, with MC-based sinogram correction, and 
with proposed geometric correction method. 
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edge artifacts, and the proposed geometric correction method significantly removes 
both the blur and edge artifacts at the same time. 

For quantitative comparison, we took CT slice images reconstructed from projec- 
tions without geometric transformation as references, and employed mean Structural 
SIMilarity (SSIM) index [19] and Mutual Information(MI) [20] to assess the image 
quality of above results by measuring the similarity between them and the corre- 
sponding references. For an image X with respect to its reference Y, SSIM index is 
defined as follows 



where x and y refer to a local window in the image X and Y respectively, \i x (fi y ) is 
the mean while a x (a y ) is the standard deviation over the window x (y), a xy is the co- 
variance between x and y, Ci and C 2 are two small positive constants. And the mean 
of SSIM index (mSSIM) is the average over all local windows. Mutual information is 
defined as 



where H(X) and H(Y) are the Shannon entropy of the image X and the reference Y re- 
spectively, and H(X, Y) is the joint Shannon entropy of the image X and its reference Y. 
We calculated mSSIM and MI values for the above slice images and listed the results 
in Tables 1 and 2 separately. Comparing the results, we can conclude that CT slice im- 
ages reconstructed from projections pre-processed by the proposed method achieve 
much higher quality than popular MC-based scheme. 

In-line XPCI data correction experiment 

Real in-line XPCI data was acquired from a mouse lung at X-ray imaging and bio- 
medical application beamline (BL13W1) of Shanghai Synchrotron Radiation Facility 
(SSRF). The experiment was approved by the Animal Experiments and Experimental 
Animal Welfare Committee of Capital Medical University (Beijing, China) and the 
approval ID is AEEI-2014-049. To prevent from deformation during imaging ap- 
proach, the mouse lung was fixed in 10% formalin solution, and dried in advance. 
And then it was put into a small tube, which was rolled by using Kapton film 
(Dupont, DE, USA), and placed on the rotary stage as illustrated in Figure 1. A Si 
(111) double-crystal was used to monochromatize the synchrotron X-ray beam. X- 
rays with photon energy of 18 keV were chosen to provide both good phase and ab- 
sorption contrast from the mouse lung. An X-ray sensitive CCD detector with high 
spatial resolution of 9 urn was placed 1.2 m away from the rotary stage along X-ray 

Table 1 Comparisons of reconstruction accuracies with mSSIM 

Slice No. NC MC GC~ 

1 057 069 092 

2 0.55 0.76 095 

3 0.54 0.83 096 

4 051 0.68 0.91 




(12) 



MI(X, Y) = H(X) +H(Y) - H(X, Y) 



(13) 
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Table 2 Comparisons of reconstruction accuracies with Ml 



Slice No. 


NC 


MC 


GC 


1 


1.52 


1.60 


2.15 


2 


1.63 


1.83 


2.60 


3 


1.62 


1.99 


2.73 


4 


1.50 


1.64 


2.18 



downstream to record projection images. During the CT data acquisition, the mouse 
lung was tilted around rotation axis of the stage from 0° to 180°, and 1296 projec- 
tions were collected totally with exposure time of 80 milliseconds for each one. The 
rotary stage was coarse adjusted before the imaging process. To improve the compu- 
tation accuracy, a simple background subtraction approach was employed before 
correction schemes were taken by deducting mean value of background region from 
projection images. 

Figure 6 shows one of CT slices reconstructed by FBP with MC-based sinogram 
correction method while Figure 7 is with the proposed geometric correction method. 
Six regions of interest (ROIs) are magnified to better visualize the interior details for 
comparison. Due to the benefits of correcting misplaced projection pixels by using 
the proposed method, slice images are reconstructed with much less edge artifacts 
and blurs. Compared to the slice image in Figure 6, more detailed structural features 
of alveoli, ribs and bronchus of mouse lung are better displayed in Figure 7. The sup- 
porting walls of two adjacent bronchi can be clearly identified as shown in ROI (a). 
Alveolar and bronchial walls are closed naturally and with less edge artifacts as dis- 
played in ROIs (b), (d) and (e). Structures of ribs are more sharp and without blurred 
as shown in ROIs (a), (c) and (e). And sharp and detailed image of soft tissue can be 
observed in ROI (f). 




Figure 6 Reconstructed CT slice from a mouse lung by FBP with MC-based sinogram correction 
method, (a), (b), (c), (d), (e) and (f) are magnified images of ROIs with corresponding labels in left 
slice image. 
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Figure 7 Reconstructed CT slice from a mouse lung by FBP with the proposed geometric correction 
method, (a), (b), (c), (d), (e) and (f) are magnified images of ROIs with corresponding labels in left 
slice image. 



Conclusions 

CT reconstruction of X-ray Phase Contrast Imaging enables to investigate internal 
microstructure of biological samples with high resolution. Mechanical imperfect or 
misalignment problem makes the reconstructed slice images blurred and with edge 
artifacts, and thus destroy the features of microstructures and reduce the spatial 
resolution of in-line XPCI. To restore the images from collected projections by in- 
line XPCI, a fast geometric correction method is proposed to determine the geomet- 
ric transform parameters properly. From the results of numerical experiments on 
synthetic and real datasets, we have the conclusion that the proposed method can 
significantly improve the image quality by removing both blurring and edge artifacts 
at the same time. Geometric correction method utilizes the symmetry of projections, 
and thus provides a simple and fast scheme to correct misplaced projection data. 
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